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ABSTRACT 

Aims. Most of our knowledge of the circumnuclear matter in Seyfert galaxies is based on the X-ray spectra of the brightest Compton- 
thick Seyfert 2 galaxies. The complete obscuration of the nuclear radiation in these sources allows us to study all the components 
arising from reprocessing of the primary continuum in the circumnuclear matter in detail, while they are heavily diluted in unobscured 
sources, often down to invisibility. 

Methods. We present the XMM-Newton RGS and EPIC pn spectra of a long (=; 100 ks) observation of one of the soft X-ray brightest 
Compton-thick Seyfert 2 galaxies, NGC 424. As a first step, we performed a phenomenological analysis of the data to derive the 
properties of all the spectral components. On the basis of these results, we fitted the spectra with self-consistent photoionisation 
models, produced with CLOUDY. 

Results. The high-energy part of the spectrum is dominated by a pure neutral Compton reflection component and a neutral iron Ko- 
line, together with Ka emission from neutral Ni, suggesting a significant Ni/Fe overabundance. The soft X-ray RGS spectrum comes 
mostly from line emission from H-like and He-like C, N, O, and Ne, as well as from the Fe L-shell. The presence of narrow RRC from 
O VIII, O VII, and C VI, the last two with resolved widths corresponding to temperatures around 5-10 eV, is a strong indication of a gas 
in photoionisation equilibrium, as confirmed by the prevalence of the forbidden component in the O VII triplet. Two gas phases with 
different ionisation parameters are needed to reproduce the spectrum with a self-consistent photoionisation model, any contribution 
from a gas in collisional equilibrium being no more than 10% of the total flux in the 0.35-1.55 keV band. When this self-consistent 
model is applied to the 0.5-10 keV band of the EPIC pn spectrum, a third photoionised phase is needed to account for emission lines 
with higher ionisation potential, although Ka emission from S XV and Fe XXVI remains under-predicted. 

Key words. Galaxies: active - Galaxies: Seyfert - galaxies: individual: NGC424 



- *~J ■ 1 . Introduction 1 kpc) high ionisation nuclear emission line region (HINER), 

/\ 70% of which was however unresolved (<200 pc) have been 

»H ; The X-ray spectrum of Compton-thick Seyfert 2 galaxies is dom- discovered (Murayama et aL 1998) . HST/WFPC2 data showed 

mated by reflection components, from both cold and ionised cir- ±e nce of a dust lane across the central t of the al 

cumnuclear matter (e.g. Matt 2000). The complete obscuration (MaJkan et a , 1998)> which heJ Mn ±e observed 

of the nuclear radiation, at least up to 10 keV, permits a clear Ay ^ 4 tQ the NLR (M ama et aL 1998) . The IRAS colours 

view of these components that are heavily diluted in unobscured ^ ke warm> suggesting that the IR emission is dominate d by 

sources, often down to invisibility. Most of our knowledge of the dugt essin of the nuclear radiation (e . g . Matt et al. 2000). 
circumnuclear matter, at least as far as their X-ray properties are 

concerned, is based on the brightest Compton-thick sources, like In X-rays, the ASCA spectrum shows a prominent iron line, 

Circinus (e.g. Matt et al. 1999; Sambruna et al. 2001; Molendi a flat spectrum, and a high [O Ill]/F(2-10 keV) ratio, suggesting 

et al. 2003), NGC 1068 (e.g. Kinkhabwala et al. 2002; Matt et al. that the nucleus of NGC 424 should be absorbed by Compton- 

2004) and Mrk 3 (e.g. Sako et al. 2000; Bianchi et al. 2005; thick matter (Collinge & Brandt 2000). This result has been 

Pounds & Page 2005). fully confirmed by the BeppoSAX observation (Matt et al. 2000; 

NGC 424 (a.k.a. Tololo 0109-383, z=0.01 17) is one of the Iwasawa et al. 2001), which measured an absorbing column of 

brightest Compton-thick Seyfert galaxies (e.g. Matt et al. 2000). about 2xl ° 24 cm \ for an estimated nuclear 2-10 keV lumi- 

Broad Ha and Hf3 lines were observed in polarised light (e.g. nosity of about 10 43 erg s NGC 424 was then observed by 

Moran et al. 2000), where the polarisation degree is about 4% af- Chandra and XMM-Newton (Matt et al. 2003). Both observa- 

ter correction for starligh. Fe II emission and an extended (about tion s were rather short (less than 10 ks each), but good enough 

to derive a few interesting properties of the source. The nuclear 

Send offprint requests to: Andrea Marinucci radiation was found to be reflected by both cold material (prob- 
e-mail: marinucci@f is . uniroma3 . it ably the inner wall of the torus) giving rise to the Compton 
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reflection component and the iron Ka fluorescent line, and by 
ionised matter, responsible for the soft X-ray emission. All these 
components are seen through a dust lane, responsible for the 
Balmer decrement and the absorption of the soft X-ray emission. 
Extended emission features in the Chandra observation are dis- 
cussed in Matt et al. (2003). We present here a new, long (~ 100 
ks) XMM-Newton observation of NGC 424, to investigate the 
properties of its circumnuclear regions in greater detail. 

2. Observations and data reduction 

The XMM-Newton observation analysed in this paper was per- 
formed on 2008 December 7 (obsid 0550950101), with the EPIC 
CCD cameras, the pn (Striider et al. 2001) and two MOS (Turner 
et al. 2001), operated in large window and medium filter, and the 
RGS cameras. Source extraction radii and screening for intervals 
of flaring particle background were performed with SAS 10.0.0 
(Gabriel et al. 2004) via an iterative process which leads to a 
maximization of the Signal-to-Noise Ratio (SNR), similarly to 
that described in Piconcelli et al. (2004). After this process, the 
net exposure time was of about 104 ks for the pn, adopting an ex- 
traction radius of 3 1 arcsec and patterns to 4. The background 
spectrum was extracted from source-free circular regions with a 
radius of 50 arcsec. Spectra were binned in order to over-sample 
the instrumental resolution by at least a factor of 3 and to have 
no less than 30 counts in each background-subtracted spectral 
channel. This allows the applicability of the^- 2 statistics. 

The RGS spectra were reduced following the guidelines in 
Guainazzi & Bianchi (2007). Data reduction pipeline rgsproc 
was used, coupled with the latest calibration files available. We 
choose a fixed celestial reference point for the attitude solution, 
coincident with the NED optical nucleus of NGC 424. Source 
spectra were extracted in regions of the dispersion versus cross- 
dispersion and Pulse Invariant versus cross-dispersion planes, 
corresponding to 95% of the Point Spread Function (PSF) in the 
cross-dispersion direction. Background spectra have been gen- 
erated using a sub-set of blank field observations, whose back- 
ground counts matches the level measured during each individ- 
ual RGS observation. The final net exposure times are about 123 
ksforRGSl and RGS2. 

3. Data analysis 

The adopted cosmological parameters are H = 70 km s _1 
Mpc -1 , f2 A = 0.73 and O m = 0.27 (i.e. the default ones in XSPEC 
12.5.1: Arnaud 1996). Errors correspond to the 90% confidence 
level for one interesting parameter (A^- 2 = 2.7), if not otherwise 
stated. The RGS spectra were not re-binned and were analysed 
using the Cash-statistics (Cash 1976). 

3. 1. RGS: phenomenological spectral analysis 

The soft X-ray spectrum of NGC 424 appears dominated by 
line emission, as commonly found in this class of sources (e.g. 
Guainazzi & Bianchi 2007). As a first step, we performed 
phenomenological fits on ^ 100-bin spectral segments, using 
Gaussian profiles at the redshift of the source (z=0.01 17), and a 
power law, both absorbed by the Galactic column density along 
the line of sight (N H = 1.8 x 10 20 cm' 2 : Dickey & Lockman 
1990). Since the model used to fit the continuum is not very sen- 
sitive to the photon index T, due to the very limited band width 
of each segment, it has been fixed to 1 . Emission lines from H- 
like and He-like C, N, O, and Ne, as well as from the Fe L- 
shell, are all clearly detected (see Table 1). The spectrum also 



NGC 424 - XMM-Newton RGS1 , O VII Ka Triplet 
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Fig. 1. RGS 1 spectrum and best fit model for the phenomenological 
analysis: the O VII triplet. The three transitions are labelled. 



presents radiative recombination continua (RRC) from Ovil, 
Ovill and Cvi. These features were fitted with the REDGE 
model in XSPEC. The width of the O VII RRC and C VI RRC 
are slightly resolved, allowing us to infer an electron tempera- 
ture of 6^3 eV for the former and 10^j eV for the latter. The 
same value inferred from the O VII RRC was adopted for the 
O VIII RRC, because it could not be significantly constrained in 
the fit. All the observed transition energies are consistent with 
the theoretical values, but there is an hint of a systemic blueshift 
in some of the strongest emission lines (e.g. C VI Ka and the 
forbidden component of the O VII triplet). On the other hand, no 
emission line is resolved, the tightest upper limits being those of 
the forbidden components of N VI and O VII Ka (<x < 220 and 
270 km s- 1 ). 

The analysis of the O VII triplet may provide precious infor- 
mation to determine whether the plasma is in photoionisation or 
collisional equilibrium. The three lines of the triplet are transi- 
tions between the n=2 shell and the n=l ground-state shell, i.e.: 
the resonance line r, corresponding to a 'So - 1 Pi transition, the 
intercombination line ; (actually a doublet: 'So - 3 P2,i) and the 
forbidden line / ('So - 3 S i). The ratio: 

G- f -^ (1) 
r 

is a good indicator of the predominant ionisation process. A 
weak resonance line compared to the forbidden or the intercom- 
bination lines is typical of plasmas dominated by photoionisation 
(G > 4). On the contrary, a strong resonance line is produced 
in collision-dominated plasmas, where G ~ 1 (e.g. Porquet & 
Dubau 2000). In the case of NGC 424, we can clearly infer that 
the gas is in photoionisation equilibrium, since G > 4.3 (see 
Table 1 and Fig. 1). 

The presence of sharp Ovil, Ovm and Cvi RRC also 
provide evidence in favour of photoionisation equilibrium. 
Radiative recombination is the capture of a free electron, to- 
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Table 1. Detected emission lines in the XMM-Newton RGS spectra of NGC 424. Energies are in keV units, wavelengths in A, fluxes in 10 5 ph 
citT 2 s -1 , kT in eV. Theoretical energies and wavelengths are from CHIANTI (Dere et al. 1997, 2009). The labelling for Fe XVII lines follows that 
of Brown etal. (1998). 



Line Id. 






Energy 


kT 


Fluxes 

(a) 


(b) 


(c) 


(d) 


CviKq- 


33.736 


0.367 


367 +0002 


- 




1.10 


0.5 


0.6 


NviKff 


29.534 
29.082 


0.420 (f) 
0.426 (i) 


4196 +0 0006 
0.426 


- 


< 0.8 


0.9 
0.2 


0.04 
0.03 


0.86 
0.17 




28.787 


0.431 (r) 


4305 +0 0007 

\j.-r~>\j~r _q 0003 




< 1.2 


0.3 


0.05 


0.25 


C VI RRC 


25.303 


0.490 


0.490!° : »J 


10+5 




0.8 


0.35 


0.45 


NVII Ka 


24.781 


0.500 


500 +0 005 


- 


6 +03 

"•"-0.3 


0.3 


0.21 


0.09 


VII Ka 


22.101 
21.807 


0.561 (f) 
0.569 (i) 


n sfii s+0.0002 
0.569 


- 


5 +0 - 8 
< 0.8 


3.5 
0.9 


0.4 
0.1 


3.1 
0.8 




21.602 


0.574 (r) 


5735 + ™™ 

-0.0005 




l + !H 

• -0.5 


0.81 


0.15 


0.66 


O VIII Ka 


18.968 


0.654 


0.6535_Q QQQ^ 




i r 03 

-0.3 


1.06 


0.93 


0.13 


OvilKjS 


18.627 


0.666 


0.669+°°,°! 




0.5 + »- 2 

-il.Z 


0.04 


0.02 


0.02 


FeXVIlM2 
FeXVII 3G 


17.097 
17.050 


0.725 
0.727 


o.7263+™°° 2 


_ 


0.3 + °' 2 

-0.2 


0.003 
0.02 


0.003 
0.02 


_ 


O VII RRC 
FeXVII 3F 


16.769 
16.777 


0.739 
0.739 


736+ 0003 




0-8^3 


0.6 
0.01 


0.14 
0.01 


0.47 


OvinK/3 


16.005 


0.775 


f> 77^+0.002 
u. / / J_ 002 




< 0.3 


0.04 


0.03 


0.01 


O VIII RRC 


14.228 


0.871 


871 + " ^ 

"'° -0.003 


6+ 3 * 

-3 


6 + °' 2 

-0.2 


0.59 


0.53 


0.07 


Ma TV V rv 

iNe ia js.a 


13.698 

1 ^ 1 

1 3. J J 1 


0.905 (f) 


904+0.001 

n Ql fy 




04+0.3 
n 

u.j_ 02 


0.40 
n 1 1 


0.22 
u.u / 


0.18 
u.uo 




13.447 


0.922 (r) 


n a? i +0.001 




5 +03 

o.j _ 02 


0.17 


0.1 


0.07 


NeX Ka 
Fe XVII 4C 


12.133 
12.119 


1.022 
1.023 


i 099+0.002 




7 +0 ' 2 

"• '-0.2 


0.13 


0.13 




Fe XVIII L 
FeXVII 7D 


10.524 
10.500 


1.178 
1.181 


i 177+O.OO5 
1 • 1 ' ' -0.004 




0.4!»i 


0.006 


0.006 




MgXI (f) 


9.168 


1.352 


i ^70+0.003 




0.4!" 


0.05 


0.04 


0.01 



(a) Calculated fluxes with the phenomenological analysis; (b) Line fluxes extrapolated from the RGS1+RGS2 best fit with CLOUDY; (c) Line 
fluxes arising from the component with higher photoionisation parameter (log U = 1.41^-^, logN H = 22.11+55°, ) ; W ^ me fluxes arising from 
the second photoionised phase (log U = 0.23^, log N H = 21.77^;*). 



gether with the emission of a photon with energy: 

tlLO n =E+Xn (2) 

where E is the initial energy of the electron, and Xn 1S the ioni- 
sation potential of the level into which the electron is captured. 
In the case of collisional ionisation equilibrium the ionisation 
potential Xn is comparable to the initial energy of the electron 
E =s kT. As a result, the recombination radiation is broad, with a 
width of ~ kT. In the case of photoionisation, typically kT «: Xn, 
and the recombination radiation is narrow. When observing this 
effect with the RGS on board XMM-Newton, the RRC in case 
of hot gas in collisional equilibrium is too broad to be visible. 
The presence of a narrow RRC is, therefore, a strong indicator 
of photoionisation. 



The O VII RRC (0.7393 keV) may be contaminated by the 3 
F component of the Fe XVII L emission line at 0.7390 keV, since 
the instrument resolution does not allow us to discriminate them. 
Indeed, we detect an emission line at ~ 0.7263 keV, which can be 
securely identified as a blend of the 3G and M2 components of 
the same species, Fe XVII. As already discussed in Bianchi et al. 
(2010), simulations with the APE C model in XSPEC showed that 
the ratio: 

Fe XVII 3F 

0.44 (3) 

Fe xvn (3G+M2) 

for a wide range of temperatures. In our case the observed ratio 
between the line detected around 0.739 keV and that at 0.7263 
keV is much higher, being 2.6^. We can therefore infer that 



4 



A. Marinucci et al.: The X-ray spectral signatures from the circumnuclear regions in the Compton Thick AGN NGC 424 



most of the flux observed at 0.739 keV is due to the O VII RRC 
component. 

3.2. RGS: CLOUDY self-consistent model 

The high quality of the RGS spectrum, coupled with the results 
from the phenomenological spectral analysis, encouraged us to 
build a self-consistent model able to reproduce the whole spec- 
trum, in a wavelength range from 8 A up to 35 A. We produced 
a grid model for XSPEC using CLOUDY 08.00 (last described 
by Ferland et al. 1998). It is an extension of the same model 
used in Bianchi et al. (2010). The main ingredients are: plane 
parallel geometry, with the flux of photons striking the illumi- 
nated face of the cloud given in terms of ionisation parameter 
U (Osterbrock & Ferland 2006); incident continuum modelled 
as in Korista et al. (1997) 1 ; constant electron density n e = 10 5 
cirT 3 ; elemental abundances as in Table 9 of CLOUDY documen- 
tation 2 ; grid parameters are log U = [-2.00 : 4.00], step 0.25, 
and logA^H = [19.0 : 23.5], step 0.1. Only the reflected spec- 
trum, arising from the illuminated face of the cloud, has been 
taken into account in our model. We also produced tables with 
different densities (n e = 10 3 - 10 4 cm -3 ): all the fits presented 
in this paper resulted insensitive to this parameter, as expected 
since we are always treating density regimes where line ratios 
of He-like triplets are insensitive to density (Porquet & Dubau 
2000). 

At first, we tried to fit the 8 A - 35 A spectrum using a single 
photoionised phase. The obtained fit is rather good, with most 
of the lines detected, and with a Cash/dof= 4885/4500. The 
best fit ionisation parameter is log U = O.Sl+g ?, with a col- 
umn density log Nr = 21.8^-|. Even if the fit was acceptable, 
we added another photoionised phase, to improve the residuals, 
for a final Cash/dof= 4865/4497. The photoionisation parame- 
ters for this best fit are log Ui = 1.41+°;°* and log C/ 2 = 0.23+°;°*, 
while the column densities values are logN Hl = 22.11 + ° °j and 
log N Hl = 21.77+°^. These values agree with similar studies 
performed on other Seyfert 2 objects (e.g. Kinkhabwala et al. 
2002; Schurch et al. 2004; Bianchi et al. 2010). The total flux 
in the 0.35 - 1.55 keV band is (2.0+°; 2 ) x 10- 13 erg cirr 2 s" 1 , al- 
most equally distributed between the two photoionised phases. 
A systemic blueshift of -230+ j° km s _1 is required by the fit, 
confirming the results found in the previous section. No signifi- 
cant improvement is obtained if different systemic velocities are 
considered for each of the photoionised phases. Similar veloc- 
ity blueshifts have been measured in the RGS spectrum of NGC 
1068, as discussed in Kinkhabwala et al. (2002). 
We tried to convolve our best fit model with a Gaussian smooth- 
ing (model GSMOOTH in XSPEC), to check if the emission lines 
are affected by a systemic broadening, and we get cr = 350+™ 
km s _1 , for a Cash/dof= 4858/4496. This result agrees with 
the RGS phenomenological analysis, even if no emission line is 
individually resolved (see previous section). No significant im- 
provement is obtained if we assume a different broadening for 
each of the photoionised phases. 



1 We also used the IUE UV fluxes reported by Dunn et al. (2006) 
to estimate an a ox ^ 1.2. Since this value is strongly dependent on 
the intrinsic X-ray flux of NGC 424, which is not directly observed, 
we decided to adopt a well-known standard continuum model, with an 
a ox = 1 .4, which is not significantly different from the one suggested 
by the data for this source. 

2 Hazy 1 version 08, p. 67: http://viewvc.nublado.org/ 
index .cgi/tags/release/c08.0 0/ docs/hazyl_08 . 
eps?revision=2 342&root=cloudy 



The phenomenological analysis agrees with the self- 
consistent model fits, and the inferred fluxes of all the emission 
lines reproduced by the CLOUDY model are shown in Table 1 . 
The two different components of the self-consistent model give 
different contributions to the whole spectrum. We can see in 
Table 1 how the main contribution to the O VII Ka line fluxes 
arises from the component with lower photoionisation param- 
eter. However, in the case of Ne X Ka, the predicted flux falls 
short the observed one. We tried to introduce a further pho- 
toionised phase with a higher photoionisation parameter, but the 
fit is not sensitive to the new component and residuals due to 
Nex Ka line are still present. At almost the same energy, an 
emission line from FeXVll is expected to be strong in colli- 
sionally ionised plasma. Since other Fe XVII lines are indeed de- 
tected in the spectrum, and their fluxes are all underestimated by 
our model, we tried to introduce a collisional component. The 
Cash/dof is = 4854/4494 and the resulting kT is 0.54+°^ keV. 
The FeXVll lines are marginally fitted, and the 0.35-1.55 keV 
flux due to the collisional phase is less than the 10% of the total 
flux. We conclude that there is no strong evidence for a colli- 
sional component in our spectrum. 

It is also interesting to note the underprediction by the model 
of the Ovil K/3 line and of the components of the Nvi and 
Ovil Ka triplets, suggesting that the resonant lines are satu- 
rated. These inconsistencies could be due to the absence of tur- 
bulence velocities in our model. It was shown that the optical 
depths to resonant absorption can be more than 10 times lower 
when the gas is affected by strong turbulence (greater than 300 
km s _1 : Nicastro et al. 1999). This effect would prevent the reso- 
nant lines to saturate up to higher column densities than in a gas 
not affected by turbulence, as assumed in our simple model. 

Apart for these lines, the residuals of the model with two 
photoionised phases are rather good, with the exception rep- 
resented by some positive residuals around 13 and 33 A (see 
Fig. 2). Since the latter are significant only in the RGS1 spec- 
trum we performed a separate analysis between the RGS spec- 
tra. It leads to a line detection only in the RGS1 spectrum 
(3.8 +2 ;° x 10~ 5 ph crrr 2 s" 1 at 0.3803 +0 ,; ,°° 1 4 keV), while in the 
RGS2 only an upper limit is obtained when the line energy is 
fixed at 0.3803 keV (< 0.6 x 10~ 5 ph crrr 2 s" 1 ). The same ap- 
proach has been used to treat the residuals around 13 A, we 
found in the RGS2 spectrum a line at 0.947+°°° 4 keV with a 
relative 0.3+°- 2 x 10" 5 ph cm s flux. The latter result is fully 
consistent with a defective pixel in the RGS2 camera while since 
the identification of the former emission line would also be prob- 
lematic, we conclude that the two previous detections are inse- 
cure 3 . 

3.3. EPIC pn spectral analysis 

3.3.1. 5-10 keV phenomenological spectral analysis 

Prominent emission lines can be clearly observed in the hard part 
of the spectrum (5-10 keV). Following Iwasawa et al. (2001) and 
Matt et al. (2003), it has been fitted with a model composed of a 
strongly absorbed (^ 10 24 cm' 2 ) power-law with F = 2, a pure 
cold reflection component (model pexrav in XSPEC, with the 
cosine of the inclination angle fixed to 0.45), and as many emis- 
sion lines as required, all of them described by Gaussian pro- 



3 We note here that a line wavelength does correspond 
to a known defective pixel in one of the RGS cameras 
(http : / /xmm .esac.esa.int/ external/xmm_user_ 
support /document at ion /uhb /node 5 9 . html#317 7 ) 
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NGC 424 - XMM-Newton RGS1 +RGS2 

Two Photoionised phases (logU=1.41 - logN H =22.11; logU=0.23 - logN H =21.77) 
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NGC 424 - XMM-Newton RGS1 +RGS2 

Two Photoionised phases (logU=1.41 - logN H =22.11; logU=0.23 - logN H =21.77) 



Energy (keV) 




20 25 30 

A (A) 



35 



Fig. 2. RGS1 (black) and RGS2 (red) self-consistent best fit (see text for detail). The 8-20 A spectrum is shown in top panel, the 20-35 A one in 
bottom panel. 



files. The strongest detected emission line is the neutral Fe Ka 
at 6.400 keV 4 , accompanied by a much fainter Fe KJ3 at 7.058 
keV. Ka emission from neutral Ni at 7.472 keV and Fe XXVI at 



4 Many of the following emission lines are indeed doublets, with 
a fixed intensity ratio between the two components of 1:2, and the 
given energies are the weighted mean (based on the values reported by 
Bearden 1967). As already shown by Yaqoob et al. (2001) and Bianchi 
et al. (2005), given the small energy distance between the two compo- 
nents, the width of the line is unaffected by the modelling with a single 
Gaussian with current X-ray instruments. 



6.966 keV are also detected, while only upper limits are found 
for the main components of the Fe XXV Ka triplet (see Table 2). 

While all other lines are unresolved, the neutral Fe Ka line 
width is <x = 81 + 12 eV. Before discussing the physical meaning 
of this width, we investigated any possible instrumental effects. 
We analysed the calibration observations performed before and 
after our observation (obsid 0412990601 and 0510780201), and 
modelled the data with a power-law (F = 1, any other choice 
does not affect the results) and three Gaussian lines, for the 
Mn Ka doublet (the separation of the lines fixed to be 0.0111 
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keV, and their flux ratio 1:2, Bearden 1967) and the KJ3. All 
line widths were kept fixed to 0. The whole model was con- 
volved with a Gaussian smoothing (model GSMOOTH in XSPEC), 
with free sigma. The best fit values for the GSMOOTH cr were 
62 + 6 and 48+6 eV, for the two data-sets. We therefore adopted 
GSMOOTH and the mean cr = 55 eV to correct the pn spectrum 
of NGC 424. 

With this correction, the line width of the neutral Fe Ka line 
is still resolved (cr = 48^0 e V). We then added a Compton 
Shoulder (CS) redwards of the line core, as expected on theo- 
retical ground, modelled as another Gaussian line, with energy 
fixed at 6.3 keV and cr = 40 eV (e.g. Matt 2002). Although the 
CS is not required by the data (A^- 2 = 1), its flux is 25 + 9% of 
the flux of the narrow core, consistent with the expectations, and 
the narrow core of the Fe Ka line is now unresolved (cr < 60 
eV). In Table 2 all the best fit results of the hard-X spectrum can 
be found and in Fig. 3 the best fit plot is shown, with a resulting 
X 2 /dof = 49/58. 

We checked the previous result by analysing the MOS spec- 
tra in the high energy band (5-10 keV). Lines' fluxes and energy 
centroids are fully consistent with the values presented in Table 
3. The width of the neutral Fe Ka line is cr < 50 eV and the 
narrow core of the line is unresolved, for a best fit ofx 2 = 35/36 
dof. 

If we consider the neutral Fe Ka line resolved, without any 
Compton Shoulder redwards of the line core, we get a cr = 48^o 
eV. With this result considerations on the geometry of the in- 
ner radius of the molecular torus can be drawn. If we assume 
that the torus has a Keplerian motion around the central super 
massive black hole, it is easy to show that the expected FWHM 
of a line produced in its inner walls should follow the relation 
2vftsini 1300(Ms/r) 1/2 sinz km s _1 , where the radius is ex- 
pressed in parsec, the mass is in 10 8 M units and i is the angle 
between the torus axis and the line of sight. The black hole mass 
of NGC 424 is estimated by means of stellar velocity disper- 
sion to be 6.02 x 10 7 M G (e.g. Bian & Gu 2007). This means 
that the expected FWHM for the iron line in this source is 1010 
r 112 sin i km s _1 . If the FWHM=5300^™ km s" 1 we measure 
with the EPIC pn is due to Doppler broadening 5 , we can estimate 
the inner radius of the torus: r = 0.04 !|j -™ sin i pc. For different 
choices of i of 30° or 60° we get central values of 0.01 and 0.03 
pc, respectively, for the inner radius. 

We can use the best fit results to provide estimates on the 
neutral gas parameters, such as the Fe and Ni abundance, and 
the ionisation state. The iron abundance is measured by the depth 
of the iron edge in the Compton reflection continuum, and it is 
measured with respect to the elements responsible for the pho- 
toabsorption below the edge, mainly oxygen and neon but with 
not negligible contributions from magnesium, silicon and sul- 
phur. Leaving this parameter free to vary in the best fit, we can 
measure it with high statistical precision, i.e. A Fe = l.OO^' 25 , in 
solar units (e.g Anders & Grevesse 1989). From the observed 
Ni Ka to Fe Ka line fluxes, we can instead estimate the rela- 
tive abundances of the two elements. The expected ratio ranges 
from 0.03 to 0.045, depending on the inclination angle and the 
incident power-law index, as discussed in Molendi et al. (2003). 
In our case, we measure 0.12 ± 0.05 (considering only the flux 
of the iron Ka core line, excluding the CS), significantly higher 
than the expected one, indicating a nickel-to-iron overabundance 
by a factor 2. This result agrees fully with previous studies on 



5 A proper modelling of the line profile should take into account that 
the rotating torus velocity distribution is double-peaked, but the quality 
of the data does not allow us to explore it in detail. 



Table 2. EPIC pn phenomenological fit (5-10 keV). See text for details. 
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high energy spectra performed on similar Compton thick Seyfert 
2 galaxies, such as NGC 1068 (Matt et al. 2004) and Mrk 3 
(Bianchi et al. 2005). The presence of an Fe K/3 line at 7.058 
keV is not statistically required, leading to an FeK/3/FeKa (core 
only) ratio < 0.10. The expected value is higher, from 0.155 to 
0.16 (Molendi et al. 2003), suggesting a moderate ionisation of 
the gas responsible for iron fluorescence. 

The detection of the Fe XXVI Ka emission line requires the 
presence of another gas with very different physical proper- 
ties from the one responsible for neutral fluorescence and the 
Compton reflection continuum. Given the limited statistics and 
the upper limits on the FeXXV Ka triplet, both a pure pho- 
toionised Compton-thin material (log U — and log/Vn = 

20.0!q3, adopting the same model described in Sect. 3.2) and 
a plasma in collisional equilibrium (kT^lO.9!^ keV, model 
APEC) reproduce the emission line equally well, with a final x 2 
statistically equivalent to the phenomenological best fit. 

3.3.2. Broad band EPIC pn spectrum (0.5-10 keV) 

Taking into account the results from the soft X-ray RGS spec- 
trum, we finally performed a self-consistent fit in the whole 
EPIC pn band (0.5-10 keV). The model for the high energy 
spectrum is the same that we adopted in the previous section, 
including a photoionisation component needed to reproduce for 
the Fe XXVI Ka emission line. As for the soft X-ray emission 
spectrum, we adopted the same model which successfully repro- 
duce the RGS spectrum, consisting in two photoionisation com- 
ponents (we fixed for all the components the same velocity shift 
measured in the RGS analysis). The overall fit is marginally ac- 
ceptable (x 2 = 256/167 dof), due to some residuals suggesting 
the presence of emission lines not reproduced by our model. 

Indeed, the addition of a line at 1.16 ± 0.01 keV is signif- 
icantly required (Ax 2 = 51), its flux of 3.4 + 0.8 x 10~ 6 ph 
ctrT 2 s _1 , agrees with the FeXVlll and FeXVll L blend ob- 
served in the phenomenological analysis of the RGS spectrum 
(see Table 1). Another emission line, at 2.44+0.02 keV, with flux 
1.1 ±0.4x 10~ 6 ph ctrT 2 s -1 , improves the fit (A^- 2 = 5). It can be 
identified with a blend of the SXV triplet (2.430, 2.448, 2.461 
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NGC 424 - XMM-Newton EPIC Pn 
High energy Best Fit 



NGC 424 - XMM-Newton EPIC MOS 
High energy Best Fit 
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Fig. 3. Best-Fit of the hard-X (5-10 keV) EPIC pn (left) and MOS (right) spectrum. 



NGC 424: XMM-Newton EPIC Pn Best fit 
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-0.01 
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Fig. 4. XMM-Newton EPIC pn 0.5-10 keV best fit and residuals. All 
the model components are shown: pure Compton reflection (grey), 
absorbed primary powerlaw (magenta), low photoionised component 
(light blue), medium photoionised component (blue), high photoionised 
component (green), Gaussian lines (red). See text for details. 



keV) 6 . More marginal detections are found at 3.13 + 0.05 keV 
(Ax 2 = 4), identified with the ArXVll Ka triplet (3.104, 3.125, 
and 3.140 keV), 3.88 + 0.04 keV (A*- 2 = 9), i.e. the CaXIX 
Ka triplet (3.861, 3.883, and 3.888 keV), and at 5.37 ± 0.07 
keV (Ay 2 = 4), consistent with Ka fluorescence from neutral 
Cr (5.412 keV). The emission line of another high ionisation 
potential species is under-predicted, i.e. Fe XXVI. This suggests 
that the gas component with the highest ionisation parameter is 



likely more complex, and cannot be easily modelled with a sin- 
gle phase, which is not able to reproduce all the observed emis- 
sion. 

The final fit is good (x 1 = 183/162 dof). The best fit values 
for the three photoionisation components are 7 : log U\ = 2.1*2'}, 
logNm > 23.3; log U 2 = 1.6+°;J, log Ate = 20.4 + 0.5; 
logt/ 3 = -1.21+°* logA/ H3 = 22.4+° 3. The i nterm ediate com- 
ponent agrees fully with the one found in the RGS analysis (see 
Sect. 3.2), both the least and the most ionised component have 
lower ionisation parameters than those extracted from the analy- 
sis in the bands where they dominate. This is due to the not neg- 
ligible contribution of the high ionisation component to the soft 
X-ray emission, reflecting the possibility that, as already men- 
tioned above, the modelling in three separate phases is too sim- 
ple, and more complex geometrical and physical dependencies 
between these components are needed in order to reproduce at 
the same time all the observed features in the broad band. 

As a self-consistent check, we extrapolated our X-ray best 
fit for the photoionised components to the optical band, in or- 
der to compare the predicted [O III] flux with the observed one. 
Our model produces an excess in this emission line, by a fac- 
tor of 16 ± 8, taking into account the errors on our best fit. The 
comparison was made with the reddening-corrected [O III] flux 
(Murayama et al. 1998), so the discrepancy is intrinsic, and not 
due to the reddening of the NLR. Moreover, no additional neu- 
tral absorption at the redshift of the source is required in our best 
fit (N H < 1.6 x 10 20 cm -2 ), at odds with the measured Ay ~1.4 
to the NLR (Murayama et al. 1998), and differently from what 
found by Matt et al. (2003), although with a much simpler mod- 
elling of the soft X-ray emission. The only contribution to the 
[O III] flux in our best fit model comes from the photoionised 
component with the lower ionisation parameter. In the X-ray 
band, this phase produces mostly O VII emission, but also strong 
emission lines from quasi-neutral S, Si, and Mg (see also Fig. 4). 



6 We excluded that this line could be an artifact of inadequate EPIC 
calibration near the gold M-edge by applying a gain fit to the data. 



7 Lower limits indicate that the parameters cannot be constrained 
above the maximum values considered in our grid model. 
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A significant contribution to these lines likely comes from the 
Compton-thick material producing the neutral iron Ka line, but 
they are completely fitted by the low photoionised phase in our 
model. Consequently, its flux may be higher than expected, pro- 
ducing part of the observed discrepancy in the [O III] emission 
line. It should also be noted that, in order to perform a detailed 
comparison between the optical and the X-ray emission of the 
NLR, a much better modelling of the spectral energy distribu- 
tion of NGC 424 should be adopted in the CLOUDY models, but 
this is beyond the scopes of this paper. 

Finally, the intrinsic 2-10 keV luminosity of NGC 424 is 
4.2 x 10 42 erg s _1 , absorbed by a neutral column density of 
1.1 ± 0.2 x 10 24 cirT 2 , it roughly agrees with the value mea- 
sured by Iwasawa et al. (2001). The 2-10 keV luminosity of the 
Compton reflection component is 9.6 x 10 40 erg s _1 , about an 
half of the reflection expected from a slab subtending at 27r. Note 
that the absorption includes Compton scattering only along the 
line-of-sight (model CABS), so the derived intrinsic luminosity 
can be considered an upper limit to the real one, depending on 
the covering factor of the absorber, and, thus, the reflection frac- 
tion is likely higher than ^ 0.5. Measuring the reflection fraction 
in Compton thick AGN is obviously impossible besides the few 
cases (such as NGC 424) in which the intrinsic continuum is 
not totally obscured below 10 keV. However 0.5 is in line 
with the typical reflection fraction of the distant reflection com- 
ponent observed in unobscured Seyfert 1 galaxies (Nandra et al. 
2007, e.g.) supporting the standard view that distant reflection in 
both obscured and unobscured AGN arises in one and the same 
medium with similar global covering factor. 

Adopting the X-ray bolometric correction presented in 
Marconi et al. (2004), we can estimate the bolometric luminos- 
ity of NGC 424 from the intrinsic 2-10 keV as L bol = 5.9 x 10 43 
erg s . This value agrees with the one that we can derive from 
the [O III] luminosity, once corrected for reddening, extracted 
from the optical data in Murayama et al. (1998), and adopting 
the bolometric correction proposed by Lamastra et al. (2009): 
Lboi — 5.1 x 10 43 erg s _1 . Assuming the BH mass already men- 
tioned in Sect. 3.3.1, the accretion rate of the source is therefore 
m ^ 6.7 - 7.5 x 10~ 3 in Eddington units. 

4. Conclusions 

We presented the XMM-Newton RGS and EPIC pn spectra of a 
long observation of the Seyfert 2 galaxy, NGC 424. The high- 
energy part of the spectrum confirmed its nature as a Compton- 
thick source, dominated by a pure neutral Compton reflection 
component, and a neutral iron Ka line with a very large EW, 
together with strong Ka emission from neutral Ni, suggesting 
a significant Ni/Fe overabundance. From the width of the neu- 
tral iron Ka line we estimated the inner radius of the torus: 
r — O.O4+0 Q2 sin i pc, where ;' is the inclination angle between 
the torus axis and the line of sight. The presence of Fe XXVI Ka 
emission can be reproduced equally well by a plasma in colli- 
sional or photoionisation equilibrium, given that only upper lim- 
its can be recovered on emission from Fe XXV. 

The soft X-ray RGS spectrum of NGC 424 comes mostly 
from line emission from H-like and He-like C, N, O, and Ne, as 
well as from the Fe L-shell, as commonly found in obscured 
AGN (e.g. Guainazzi & Bianchi 2007). The presence of nar- 
row RRC from O VIII, O VII, and C VI (the latter two with re- 
solved widths corresponding to temperatures around 5-10 eV) is 
a strong indication of a low-temperature gas, which must there- 
fore be in photoionisation equilibrium to have such an high ion- 
isation state. This is confirmed by the prevalence of the forbid- 



den component in the O VII triplet. Indeed, a self-consistent pho- 
toionisation model well reproduces the RGS spectrum, adopting 
two gas phases with different ionisation parameters. This gas is 
likely to be identified with the ionisation cones of the NLR, as 
suggested for most Seyfert 2 galaxies (e.g. Kinkhabwala et al. 
2002; Schurch et al. 2004; Bianchi et al. 2010), although a con- 
tribution from inner regions (i.e. the torus, the BLR) cannot be 
ruled out. Any contribution from a gas in collisional equilibrium 
should not exceed 10% of the total flux in the 0.35-1.55 keV 
band. 

When this self-consistent model is applied to the 0.5-10 keV 
band of the EPIC pn spectrum, a third photoionised phase is 
needed to account for emission lines with higher ionisation po- 
tential, although Ka emission from S XV and Fe XXVI remains 
under-predicted. 
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